The Journal of Physical Chemistry B

Protein translation is a highly regulated process influenced by multiple factors at the initiation, elongation, and termination stages. One notable regulatory element of the ribosome is the CAR interaction surface, a three-residue motif in the structure of the ribosome composed of C1274 and A1427 of S. cerevisiae 18S rRNA (corresponding to C1054 and A1196 in E. coli 16S rRNA) and R146 of ribosomal protein Rps3. CAR is highly conserved and positioned adjacent to the amino-acyl (A site) decoding center. It establishes hydrogen bonds with the +1 codon next in line to enter the ribosome A site, acting as an extension of the tRNA anticodon and forming base-stacking interactions with nucleotide 34 of the tRNA. However, despite CARs enzymatically strategic positioning within the ribosome, its functional relationship with the A site remains poorly characterized. Using molecular dynamics (MD) simulations, we examined the interplay between the A site and CAR site, revealing sequence-dependent modulation of H-bonding and {pi}-stacking interactions within and between the two sites. These findings highlight the interplay between the A site and CAR site, suggesting a structural and functional connection between these two regions of the ribosome that may contribute to mRNA sequence-specific tuning of translation elongation. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=91 SRC="FIGDIR/small/714784v1_ufig1.gif" ALT="Figure 1"> View larger version (22K): org.highwire.dtl.DTLVardef@1919efaorg.highwire.dtl.DTLVardef@15c4882org.highwire.dtl.DTLVardef@19c7782org.highwire.dtl.DTLVardef@16a1246_HPS_FORMAT_FIGEXP M_FIG C_FIG

Intramolecular diffusion is an important, but often overlooked, property of intrinsically disordered proteins and plays an important role in folding, assembly and aggregation. Fluorescence resonance energy transfer (FRET) is used to observe reconfiguration over nanometer length scales while close range quenching over Angstrom length scales provides a complimentary view with different dynamics. There are several probe/quencher pairs that have been employed with varying levels of quantification of the quenching rate. Here we measure the electron transfer quenching parameters of the fluorophore Atto-655 by tryptophan using fluorescence correlation spectroscopy. Measurements with varying concentrations of quencher with low diffusion yield a distance dependent quenching rate. These parameters provide for a more quantitative analysis of measurements of intramolecular diffusion, particularly in crowded environments.

Invertebrate vision relies on bistable visual pigments flipping upon photon absorption between rhodopsin and metarhodopsin states. In living butterflies, the UV-VIS absorption spectra of rhodopsin and metarhodopsin, respectively with 11-cis and all-trans isomers of 3-hydroxy-retinal (A3) chromophore, can be conveniently recorded from the eyeshine, the light reflected from the compound eye after passing twice through the light-guiding rhabdoms. * Here, a microscope coupled with a broadband LED source and a microspectrometer was used to record photorelaxations reported in eyeshine reflection spectra. Fitting temporal exponential relaxations to log-reflectance arrays yielded transient and baseline spectra that are analogous to absorbance difference and sum, respectively. Both types of spectra were subjected to singular value decomposition and to fitting of templated visual pigment absorption spectra. * The compound eye of the high brown fritillary Fabriciana adippe was exposed to a series of second-long broadband light pulses, causing photorelaxations with time constants between 40 and 120 ms that led to 80% metarhodopsin in equilibrium. The transient and baseline spectra were fitted with pigment templates, estimating the alpha peak wavelength 547-552 nm for rhodopsin and 496-501 nm for metarhodopsin. The metarhodopsin to rhodopsin alpha peak absorbance ratio 1.25-1.35 is consistent with the isosbestic wavelength at 530 nm. The second isosbestic wavelength indicates that rhodopsin beta (UV) peak absorbs more strongly than metarhodopsin below 405 nm. * Baseline spectra, which were not explicitly analysed in previous studies, enable concatenation of exposures, monitor long-term changes of pigment, and enhance the estimation of beta peak parameters. * The method can be directly used in many butterflies and could be adapted to other insects, particularly fruitflies, facilitating studies of the relation between the visual pigment spectra and the opsin sequences. Spectroscopic results can be complemented with physiologically measured photoreceptor spectral sensitivity datasets and analysed with the same global fitting procedure.

The conformational dynamics of a model cationic protein in water and in the presence of anionic sodium dodecyl sulphate (SDS) and cationic cetyltrimethylamonium bromide (CTAB) surfactants at different concentrations were investigated using all-atom molecular dynamics simulations. Free-energy landscapes constructed along principal components reveal a compact, well-defined native basin at 25 {degrees}C in water, whereas elevated temperature (100 {degrees}C) induces a broadening of the conformational space and the emergence of multiple metastable states. The presence of surfactants further modulates this behavior in a concentration-dependent manner. Cluster population analysis shows that SDS promotes a highly heterogeneous ensemble characterized by reduced dominance of the native-like cluster, while CTAB partially protects the protein from thermal denaturation at higher concentrations. Radial distribution functions demonstrate strong accumulation of SDS headgroups around the protein and pronounced insertion of SDS alkyl tails into hydrophobic protein regions, indicating direct hydrophobic destabilization and micelle-assisted unfolding. In contrast, CTAB exhibits weaker headgroup association owing to electrostatic repulsion and reduced tail-hydrophobic contacts, suggesting a less disruptive interaction mechanism. At high concentration, CTAB aggregates provide a structured hydrophobic environment that stabilizes the folded state and suppresses denaturation. Together, these results provide a molecular-level picture of how surfactant chemistry and concentration govern the conformational stability of a cationic protein, highlighting the dominant role of hydrophobic interactions in surfactant-induced denaturation at high temperature. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/717321v1_ufig1.gif" ALT="Figure 1"> View larger version (24K): org.highwire.dtl.DTLVardef@dcf96aorg.highwire.dtl.DTLVardef@17acdc7org.highwire.dtl.DTLVardef@15bdc2borg.highwire.dtl.DTLVardef@1d39f3c_HPS_FORMAT_FIGEXP M_FIG C_FIG

The molten globule (MG) state is an intermediate in the unfolding pathway of proteins, typically triggered by denaturing agents such as urea, extreme pH, high pressure, or heat. The microscopic details of such states are far from understood. Here we study the MG states in protein Hen Egg-White Lysozyme (PDB ID: 1AKI) using microscopic constant pH molecular dynamics (CpHMD) simulations and experiments across a wide pH range. We observe that the titratable residues act as key drivers of conformational fluctuations, promoting the emergence of MG states at extreme pH. These states display partial unfolding, and small global structural changes (< 7% deviation). Hydration around the fluctuating acidic residues shows reduced water density and weakened hydrogen bonding at low pH. At high pH, hydration around acidic residues increases relative to pH = 7, whereas hydration around basic residues decreases. The translational and rotational dynamics of hydration water also exhibit pronounced pH dependence: the translational diffusion coefficient (Dtrans) increases linearly with decrease in pH in acidic medium and increases linearly with increasing pH in the basic regime. The rotational diffusion (Drot) shows similar dependencies on pH except a break at pH {approx} 4 corresponding to acidic residue pKa values. Our results may be useful to identify ligand binding of lysozyme in extreme pH conditions.

Single molecule methods have become prevalent tools in elucidating molecular processes across various life science fields over the past three decades, driving breakthroughs in understanding their underlying molecular mechanisms. In our study, we employed two single-molecule methods, Forster Resonance Energy Transfer (smFRET) and high-resolution optical tweezers, to investigate the transcription of Mycobacterium tuberculosis RNA polymerase (MtbRNAP) from initiation through to termination. We aim to provide a set of comprehensive biophysical tools to deepen our current understanding of MtbRNAP and its transcription factors. These experimental assays represent an important step towards unraveling the molecular dynamics and interactions that support transcription in Mycobacterium tuberculosis.

Amyloid beta (A{beta}), one of the hallmark proteins of Alzheimers Disease (AD), aggregates into plaques that are strongly linked to cognitive decline and neuronal death. Reducing its aggregation propensity may provide a strategy to slow the progression of AD. While chirality modulation has emerged as an innovative approach to disrupt this process, research has primarily focused on alterations at the C position, often overlooking the impact of the second chiral center, such as the C{beta} atom of Threonine. Furthermore, the underlying mechanisms governing these chiral effects remain elusive. Given the intrinsically disordered nature of the A{beta} peptide, we employed temperature-replica exchange molecular dynamics (T-REMD) simulations to explore its rugged conformational landscape. We considered sequence mutations (A2T, A2V), N-terminal chirality inversion of the first six residues (A2V1-6D and WT1-6D), and alteration of the second chiral center (C{beta}) of Threonine (A2TC{beta}). By analyzing the effect size and population change induced by these mutations and chiral modulation, we concluded that the modulation at the N-termini is not confined locally but also exerts specific effects on the central hydrophobic core (CHC) region. Inspection of their free energy landscape and representative structures reveals that the protective or pathogenic effects of these variants correlate with their similarity to the wild type (WT) ensemble. Beyond these static thermodynamics analyses, a direct connection to phase transitions was made by estimating heat capacity as a function of temperature. Both analyses predict that A2TC{beta} may exert a pathogenic effect, in contrast to the protective nature of A2T. These findings offer a deeper understanding of the effects of site-specific mutations and chirality and shed light on the development of advanced therapeutic strategies for AD.

We have conducted all atom molecular dynamics simulations of POPC and DPPC lipid bilayers using AMBER Lipid21 force field with eight different water models, including SPC/E, TIP3P, TIP3P-FB, TIP4P-FB, TIP4P-Ew, TIP4P/2005, TIP4P-D, and OPC, to identify the most compatible one without any modification. A number of parameters have been computed in order to understand the structure of the lipid bilayer: Area per lipid, Isothermal compressibility modulus, average Volume per lipid, electron density profile, bilayer thickness, X-ray and neutron scattering form factors, deuterium order parameter, and radial distribution function. The estimated Area per lipid, Isothermal compressibility factor, volume per lipid and bilayer thickness are highly consistent with experimental results for the SPC/E water model, indicating its suitability with the AMBER Lipid21 force field, insted of any modification. The bilayer electron density profiles of both the lipid bilayers demonstrate a little augmentation of water penetration with respect to the membrane surface for TIP4P-D water model. However, the experimental X-ray and neutron scattering form factors are aligning well with the simulated results for all studied water models, and TIP4P-D shows better for X-ray data. The deuterium order parameter for lipid acyl chains value less than 0.25 for all observed water models, depicting their disorderness for both the lipid bilayers. The lateral diffusion and reorientation autocorrelation function of the lipid molecules in both the bilayers are computed to reveal their dynamics across all water models. In comparison to other water models, the simulated trajectories predict better structure and reasonably fair dynamic properties for the SPC/E water model. The TIP4P-Ew water model reproduces the lateral diffusion co-efficient in close agreement with experiment. Reorientational dynamics for both the lipids in the bilayers for eight different water models are observed; the presence of slow and slowest time components corresponds to the lipid axial motion (wobble motion) and Twist/Splay motions. So, in view of the overall performance of the different water models with the AMBER Lipid21 all atom force field in reproducing membrane physical properties, the SPC/E water model appears to be an optimal choice.

Mycobacteria are responsible for causing severe illnesses like tuberculosis and leprosy in humans. Studying the mycobacteria cell envelope presents a significant challenge due to its intricate lipid compositions and structural variations and also its harmful nature in a typical experiment setting. In this study, we use all-atom molecular dynamics simulation to study mycobacterial inner membranes (MIMs). By incorporating different types of phosphatidyl-myo-inositol-mannosides (PIMs) and their glycoconjugates such as lipomannans (LM) and lipoarabinomannans (LAM) lipoglycans, we have constructed both symmetric and asymmetric membrane systems to study the MIM structure and dynamics under varying compositions of each lipid type. Our results show that the phospholipid/PIM-rich inner leaflet remains a stable, fluid bilayer, and the outer leaflet structure and dynamics are heavily governed by lipoglycan surface density. Importantly, as LM/LAM concentration increases, the polysaccharide chains shift from flexible, membrane-lying orientations to a compact brush-like state aligned with the membrane normal. This crowding significantly reduces the solvent-accessible volume and limits direct interactions between LM/LAM sugars and the outer leaflet surface. Furthermore, we observe that high lipoglycan presence in the outer leaflet slows lipid diffusion across the entire bilayer, demonstrating a dynamic coupling between the two leaflets. By resolving these LM/LAM sugar-level dynamics and their impact on membrane-wide properties, this study provides a molecular framework for future MIM modeling and simulation with various (peripheral) membrane proteins to better understand how the MIM functions as a regulated physical barrier and a platform for mycobacterial virulence.

The catalytic mechanism of the hairpin ribozyme has remained controversial for more than two decades, with different experimental approaches often supporting distinct mechanistic interpretations. In this work, we investigate the conformational landscape of the active site along several proposed reaction pathways using all-atom molecular dynamics simulations in explicit solvent combined with enhanced sampling techniques. Specifically, we employ Hamiltonian replica exchange simulations to extensively explore active-site conformations without relying on predefined collective variables, enabling a broad characterization of the structural ensembles associated with multiple protonation states along three candidate reaction mechanisms. Our simulations suggest that a dianionic general acid/general base pathway involving direct participation of A38 and G8 is unlikely to proceed through well-defined intermediates with catalytically competent geometries. In particular, states associated with G8 deprotonation and subsequent O2 deprotonation exhibit strongly distorted active-site arrangements that appear poorly suited for reaction progression. Although highly synchronous protontransfer steps cannot be excluded, the required deprotonation of G8 remains difficult to reconcile with neutral pH conditions. In contrast, monoanionic pathways in which the non-bridging oxygens of the scissile phosphate act as transient proton relays produce intermediates that sample geometries favorable for the nucleophilic addition and leaving-group elimination steps of the reaction. These mechanisms do not require direct catalytic involvement of G8 while remaining compatible with potential acid catalysis by protonated A38+. Our results provide a unified conformational perspective on competing mechanistic scenarios. The ensembles generated here offer a foundation for future QM/MM and ML/MM calculations aimed at quantitatively resolving the free-energy landscapes governing hairpin ribozyme catalysis. Finally, the present strategy could easily be applied to other biomolecular systems with high conformational plasticity, including other ribozymes.

Like other proteins, monoclonal antibodies - important biodrugs- are subject to post translational modifications, especially the N-glycosylations. However, the effect of the N-glycosylations remains poorly studied and atomistic details about their influence are rarely available. Moreover, the few existing studies focus on the prevalent immunoglobulin G1. To go further in the understanding of the impact of glycosylations, we have carried out a comparative exploration of the effect of N-glycosylations on two different classes of antibodies, namely Mab231, an IgG2 and the pembrolizumab, an IgG4. The two antibodies differ by their sequences, their length, their 3D structure but also by the location and composition of the glycans. In the present work, detailed and important information were gained through molecular dynamics simulations where both monoclonal antibodies were studied without and with the presence of their glycans. The results of 1.5 {micro}s of sampling for each system show that glycosylation does not drastically alter the overall conformational landscape of either antibody, whatever the metrics considered. However, it measurably modulates local flexibility, inter-domain correlated motions, and the relative orientation of the Fab arms with respect to the Fc domain, with statistically significant shifts in key geometric descriptors. Importantly, contact analysis reveals that glycan interactions extend beyond the Fc region to reach Fab residues. The allosteric network calculations demonstrate that the influence of Fc-bound glycans propagates even until the Fab framework regions in both mAbs, which could impact the antigen binding. The nature and magnitude of these effects are subclass-dependent, reflecting differences in glycan composition, hinge architecture, and three-dimensional organization Our findings challenge the prevailing view that Fc glycosylation uniformly promotes CH2 domain opening. More importantly, it underscores the necessity of considering full-length structures and IgG subclass diversity in glyco-engineering strategies.

Intrinsic tryptophan fluorescence is widely used as a sensitive reporter of protein conformational dynamics, yet the molecular origin of its temperature-dependent modulation remains unclear. Here we investigate the conformational dynamics of Trp134 in bovine serum albumin (BSA) using molecular dynamics (MD) simulations, free-energy calculations based on umbrella sampling and WHAM, quantum mechanical (QM) calculations, and QM/MM approaches. MD simulations show that the global structure of BSA remains stable while temperature induces a gradual population shift from the Ia+ to the Ia- rotamer. The corresponding free-energy landscapes reveal that this shift arises from subtle changes in basin stability and transition barriers along the rotameric coordinate. In contrast, standalone QM calculations on isolated tryptophan predict different energetic trends, highlighting the sensitivity of rotamer stability to electronic-structure treatments and environmental effects. QM/MM calculations partially reconcile these differences by incorporating the protein environment. Together, these results suggest that temperature reshapes the rotamer free-energy landscape of Trp134, leading to population shifts that modulate intrinsic tryptophan fluorescence in proteins.

There has been a surge in antibiotic resistance in recent years, making traditional antibiotics less effective against key pathogens. RNA has recently emerged as a potential target for antibiotics due to its involvement in crucial microbial functions. It is possible to expand the range of therapeutic targets by using RNA-based therapies, but it remains necessary to improve the molecular-level understanding of interactions between RNA and known and potential binders. The SAM-I riboswitch, which controls the transcriptional termination of gene expression involved in sulfur metabolism in most bacteria, is an excellent ligand target. Thus, understanding its behavior with and without ligand complexes would be very helpful for drug design applications. In this manuscript, we studied the interactions between the SAM-I riboswitch and its natural ligand, SAM, which controls riboswitch function, and compared those interactions to its interactions with the very similar small molecular SAH, which does not control riboswitch function, and to its interactions with a potential binder JS4, identified via virtual screening. From our simulations, we gain a deeper understanding of small molecule interactions with the SAM-I riboswitch. The results reveal how differently the small molecules (SAM, SAH and JS4) bind to and potentially induce conformational changes in the riboswitch. Our findings offer valuable insight into the molecular mechanisms underlying riboswitch RNA-ligand interactions for the design of more effective RNA-targeting therapeutics.

Proteins in solution adsorb to the corona of nanoparticles such as single-walled carbon nanotubes (SWCNTs), but these interactions are difficult to predict and analyze due to ambiguities in the structure of the latter. In this work, we employ ss(GT)15-DNA wrapped SWCNTs, a commonly used fluorescent sensor construct, to examine protein adsorption by quantifying binding dissociation constants and characterizing the corresponding photophysical effects. A library of 20 proteins are used to evaluate adsorption-induced changes in photoluminescence (PL) intensity ({Delta}I/I0) and emission wavelength upon solution phase binding. We find that 15 proteins produce monotonic dose-response behavior well described using a single-site Langmuir model. Alternatively, five proteins exhibited more complex, non-monotonic behavior consistent with a two-step binding model representing protein-protein interactions coupled to adsorption. The study reveals that metalloproteins, which comprised 12 of the 20 proteins in the library, induced greater PL quenching compared with metal-free proteins for this system, with maximum binding-associated quenching ({Delta}I/I0) of 94% for metalloproteins versus 20% for metal-free proteins. For metalloproteins, we introduce a proximity-based quenching framework in which protein size provides a coarse proxy for cofactor-SWCNT separation, offering a mechanistic interpretation of the observed quenching variation across proteins. Together, these results establish the use of metal coordination sites, such as those in metalloproteins, to assist the transduction of certain nanoparticle fluorescent sensors, helping with sensor probe design and interpretation in biological environments.

The repair of photo-induced DNA lesions through nucleotide excision repair machinery is still the source of important questions. It has been observed that the repair rate of the different cyclobutane pyrimidine dimers, i.e. the photoproducts induced by dimerization of two {pi}-stacked pyrimidines (T<>T, T<>C, C<>T, C<>C), depends on the nucleobases involved in the lesion. TT derivatives (T<>T) are removed more slowly than those containing cytosine, especially in 5. Using all-atom molecular dynamics simulations and free-energy calculations, we demonstrate that the variation of the repair rate observed in human skin and in cultured cutaneous cell is associated to the recognition of the four lesions by the DDB2 protein moiety, and more specifically by the differential structural deformation induced on the complementary strand. Indeed, while C<>C and C<>T induce a larger deviation on the groove parameters, T<>T and T<>C, instead, affect DNA structure to a lesser extent. less affected. These effects then hamper differentially the downstream recruitment of the repair complexes. The observed DNA deformation correlates with the experimental repair rate and provides a structural rationale for the different repair rates of CPD by nucleotide excision repair machinery. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=105 SRC="FIGDIR/small/724087v1_ufig1.gif" ALT="Figure 1"> View larger version (43K): org.highwire.dtl.DTLVardef@cf6b6dorg.highwire.dtl.DTLVardef@195e35forg.highwire.dtl.DTLVardef@1829296org.highwire.dtl.DTLVardef@165baba_HPS_FORMAT_FIGEXP M_FIG C_FIG

Protein dynamics are essential to biological function, yet understanding whether deep learning models contain information about these dynamics remains an open question. In this study, we quantitatively investigate the capacity of deep learning structure generation methods to predict protein flexibilities by directly comparing residue-level mean squared fluctuation (MSF) profiles derived from structural ensembles with experimental or simulation-informed flexibility profiles. We assembled four diverse benchmark datasets representing different types of structural information, including 70 NMR ensembles, 43 X-ray crystallographic protein pairs in two distinct conformational states, 82 high-resolution cryo-EM structures, and molecular dynamics simulations of 10 proteins. Utilizing AlphaFold3, AlphaFold2, and RosettaFold to generate multiple structural models, we applied ranksort normalization to place the profiles on a comparable scale and quantified similarity primarily using cosine and Pearson similarities. Our results demonstrate that the flexibility predictions from deep learning-generated models agree well with experimental data, suggesting that fluctuations in these predicted ensembles can serve as effective proxies for protein flexibility. Notably, AlphaFold3 consistently produced the best results across the datasets. We also observed that flexibility prediction accuracy generally improves as the number of models increases up to 15, and our findings remain robust even when terminal residues are excluded from the analysis. To facilitate broader application, we provide three publicly accessible Jupyter Notebooks to calculate MSF from deep learning outputs. Ultimately, this work provides evidence that deep learning structural ensembles can serve as proxies for protein flexibility.

Surface functionalization of inorganic quantum dot nanoparticles is of great interest in the application of these materials toward a wide range of biological applications where membrane interactions are critical. The use of amphiphilic lipids to functionalize the surfaces of quantum dots represents a promising alternative to produce water-soluble and membrane-active materials with facile tuning of the quantum dots surface properties. Here, we demonstrate an experimental approach that yields lipid-coated quantum dots with highly tunable surface charge by controlling the concentration of cationic lipids during preparation. Through fluorescence-activated cell sorting assays, we show that these cationic lipid-coated quantum dots can enhance membrane interactions and increase membrane labeling density in live HEK293 cells. We further employed coarse-grained molecular dynamics simulations to model the lipid self-assembly process using an implicit solvent force field and subsequently model the adsorption of lipid-coated quantum dots to model membranes. Our simulations show that we can control the effective surface charge of lipid-coated quantum dots and influence the strength of adsorption to oppositely charged lipid membranes, a process that is mediated by the release of counterions at the quantum dot-membrane interface. This work supports the future development of biocompatible and water-soluble inorganic nanoparticles with highly tunable surfaces, and provides mechanistic insight into how different lipids can influence nanoparticle-membrane interactions at a molecular scale.

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are are a family of voltage-gated, cyclic-nucleotide modulated Na+/K+ channels that regulate spontaneous rhythmic electrical activity in both the heart and the brain. Understanding differences in the responsiveness to cyclic adenosine monophosphate (cAMP) modulation between HCN isoforms would offer insight into the specific binding interactions that drive channel activation. Using all-atom molecular dynamics (MD) simulations and the free-energy perturbation (FEP) approach, we determined the absolute binding free energy of cAMP to the the cyclicnucleotide-binding domain (CNBD) of HCN isoforms 1-4. By studying the free-energy of ligand binding to the various isoforms of HCN, our study advances the understanding of HCN channel activation and modulation mechanisms. Overall, our work offers insight into explaining differences in channel sensitivity across the isoforms of HCN.

The apparent pKa of ionizable lipids in lipid nanoparticles (LNPs) is a key determinant of RNA encapsulation during formulation and endosomal release after cellular uptake. However, it is difficult to predict the effective pKa of a given ionizable lipid solely from its solution pKa, because it is sensitive to the membranes composition, as well as solution conditions such as the salt concentration. We developed a simple continuum electrostatics model, based on Gouy-Chapman theory, to predict the shift in effective pKa for ionizable lipids in lipid bilayers as a function of salt concentration and membrane composition. We derive equations for the surface potential and fraction of lipids charged, which are solved self-consistently as a function of solution pH to extract the titration curve and effective pKa. The model shows that the shift in effective pKa is largest when the concentration of titratable lipid is high, and the effect is diminished by increasing salt concentration. We provide a python implementation of the model and an interactive notebook that will allow users to further easily explore the predicted pKa shifts as a function of formulation variables.

Intercalation of small molecules between DNA base pairs affects DNA conformation, disrupting essential cellular processes including replication, transcription, and repair. We investigated conformational changes in 18-mer DNA upon intercalation of doxorubicin, SYBR Gold and YOYO-1 using extensive MD simulations. Two main patterns for the intercalation were identified: RISE-type intercalation occurs between adjacent base pairs and extends the DNA helix with decreased twist angles, while OPEN-type intercalation proceeds through base-pair opening without significant DNA extension. Kinetic analysis revealed that association rates for intercalation followed the order: first YO-moiety (mono-intercalation) > SYBR Gold > doxorubicin > YOYO-1 (bis-intercalation). Free energy landscape showed that forces at DNA termini reached up to 117 pN during stretching. Notably, base pairs adjacent to intercalators were protected from strand separation, accompanied by additional helical unwinding. MM-PBSA/GBSA analysis revealed that the driving force for intercalation is the stacking energy, and the binding affinity was highest for minor groove binding. Persistence length decreased with single molecule binding but recovered with two molecules due to their electrostatic repulsion. Mechanical properties of intercalated DNA showed position-dependence, demonstrating that multiple intercalation modes coexist in solution. The heterogeneous nature of intercalation explains why experimental measurements reflect ensemble averages rather than single binding configurations.